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Abstract 

The matrix element of a bound electron interacting with the nucleus through exchange of a Z 
boson is studied for the gauge invariant case of 2s±/2 ~ %Pi/2 transitions in hydrogenic ions. The 
QED radiative correction to the matrix element, which is —a/2n in lowest order, is calculated to 
all orders in Za using exact propagators. Previous calculations of the first-order binding correction 
are confirmed both analytically and by fitting the exact function at low Z. Consequences for the 
interpretation of parity nonconservation in cesium are discussed. 
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I. INTRODUCTION 



The calculation of radiative corrections in atoms with low nuclear charge Z is facilitated 
by the fact that binding corrections, which enter as powers and logarithms of Za, are rela- 
tively small, and can be treated in perturbation theory. For atoms with high nuclear charge 
the perturbation expansion converges more slowly, and for highly-charged ions the expan- 
sion is generally avoided, which is possible when numerical methods are used to represent 
the electron propagator. This approach, first introduced by Wichmann and Kroll [1] for the 
vacuum polarization and Brown and Mayers [2] for the self-energy, has been applied to the 
calculation of both energy levels, notably by Mohr and collaborators [3] , and more recently 
to matrix elements, specifically hyperfine splitting (hfs) and the Zeeman effect [4, 5]. 

It is of interest to further extend this kind of radiative correction calculation to the 
parity nonconserving (PNC) process 6si/2 — > ls\/2 in neutral cesium [6]. Corrections to 
this transition are of importance for the question of whether a breakdown of the standard 
model is present for cesium PNC. Specifically, if the radiative correction to the electron-Z 
vertex is taken to be its lowest-order value, — a/27r, then based on the the present status of 
other corrections to PNC which have included a number of significant shifts only recently 
considered that arise from the Breit interaction [7] and vacuum polarization [8], a discrepancy 
with experiment of approximately la would result. Given the presence of other indications 
of possible problems with electroweak tests of the standard model, specifically the NuTev 
result [9] and hadronic asymmetries in Z — > bb [10], a discrepancy in cesium PNC could be 
an indication of new physics. 

However, it is known that binding corrections to the similar matrix element involved in 
hfs are very large for highly- charged ions. That this is so is not surprising, given the first 
two terms of the one- loop vertex correction to hfs [5], 



where Ep is the lowest-order hfs energy. Already at Z = 9 the leading binding correction 
leads to a change in sign of the hfs, and at Z = 55 the formula would predict —2.72 - Ep, 
as compared to the low-order, uncorrected value of +0.5^ Ep. Of course, with Za = 0.4, 
the above equation, even with known higher-order terms included, cannot replace an exact 
evaluation. As mentioned above, such evaluations have been carried out by a number of 
groups, and the complete answer turns out to be —3.02 ^ Ep [5]. 




(1) 
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It is possible to carry out a parallel analysis for radiative corrections to PNC. If we define 
the lowest-order PNC matrix element as Qo and the one-loop radiatively corrected matrix 
element as Qr, with 

Q R = -Q R(Za), (2) 

71 

the first two terms of R(Za) are 

R(Za) = ~-(2hi2 + ^yZa, (3) 

where the first term is part of the standard radiative correction for atomic PNC [11] and 
the leading binding correction was first calculated in Ref. [12]. For the case of cesium this 
formula changes the coefficient of a/n from -0.5 to -2.98, changing a negligible -0.12 percent 
to a significant -0.69 percent shift. This largely removes the 2a discrepancy between theory 
and experiment. 

There are a number of issues that must be addressed before accepting the -0.69 percent 
shift at face value. Firstly, just as with hfs, an approach that does not rely on expansion 
in Za is required. Even though the first two terms in Eq. (1) for the vertex correction to 
hfs give an answer within 12 percent of the total answer, there is no reason we know of for 
this to be true in general. Secondly, it is not clear that it is correct to use Z = 55 in the 
above equation. When the cesium 6s Lamb shift, which is also governed by short distance 
effects, is studied with all-orders methods [13, 14], a much smaller effective nuclear charge 
is seen, specifically about 14. Thirdly, an important difference between PNC and hfs is the 
role of gauge invariance. In the latter case the initial and final states are real physical states. 
However, the Z boson vertex does not involve two physical states, instead involving either a 
6si/ 2 or 7s 1/2 state and an intermediate state with py 2 quantum numbers. While it can be 
shown that Eq. (3) is still valid in this case, higher-order binding corrections will be gauge 
dependent. 

To address the last issue, we choose here to work with a gauge-invariant quantity, the 
matrix element of the weak Hamiltonian 

G 

H w = Qw-j=lolsPN{r) (4) 

between the 2si/ 2 and 2p ± / 2 states of a hydrogenic ion, where pjv(r) describes the distribu- 
tion of the weak nuclear charge, which is close to the neutron distribution. While a finite 
distribution will be used for pAr(r), the atomic 2si/ 2 and 2pi/ 2 states will be chosen to be 
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solutions of the Dirac equation with a point nucleus, so the energies of these two states are 
equal. This allows radiative corrections to PNC to be studied nonperturbatively to all orders 
in Za in a manner parallel to that used for hfs [5, 15], and in particular gives information 
about the Za behavior of the function R(Za) that will be useful when the cesium problem 
is addressed, as will be discussed in the conclusion. 

The plan of the paper is the following. The lowest-order matrix element Qo is treated 
in Sec. II. In Sec. Ill we give a derivation of the radiative correction formulas, and in 
Sec. IV evaluate R(Za) to first order in Za, confirming the result of Ref. [12]. In Sec. V we 
rearrange the formulas in a way that allows for an exact numerical evaluation, and present 
the details of such a calculation for the range Z = 10 — 100. In the last section, it is shown 
that the numerical evaluation at low Z agrees with the perturbative expansion, and the 
higher-order binding corrections inferred. Prospects for extension of the calculation to the 
actual experiment, where a laser photon is present driving the 6si/ 2 — 7s 1 / 2 transition, are 
also discussed. 



II. LOWEST-ORDER CALCULATION 



The matrix element of the weak charge operator in lowest order is 

Qo = Qwv = / ci 3 r^(r)7 5 ^(f)pAr(f), (5) 



where we shall from now on suppress the overall factor QwGfI 'a/8, use w to denote the 
2pi/2 state, and v the 2si/ 2 state. The nuclear distribution is chosen to be uniform, with a 
radius i?o fixed so that the root-mean-square radius agrees with a fermi distribution with a 
thickness parameter 2.3 fm and a c parameter given in Table I. Because of the simplicity of 
the uniform distribution the matrix element can be evaluated analytically, and is 

QiZ 4 a f2ZR \ 2 ~<- 2 ^1 + 2^ _2Mo , 



Here 7 = y 1 — (Za) 2 , N 2 = y 2(1 + 7), and a is the Bohr radius. We note the singularity 
of this expression as Ro — > 0, which at small Z manifests itself as a logarithmic dependence 
on R , as can be seen from the Taylor expansion in Za of the above, 



^ y/3iZ 4 a r r lnr N2 / . 55 x\ „,„ s 

Qo = 32^ 6 i 1 + {Za) \ lnX ~ lE + 24 " 8 J + ( } 



(7) 



where x = ZR /a and 7^ = 0.577. . . is Euler's constant. Results of Q are tabulated in 
Table I. 



III. DERIVATION OF RADIATIVE CORRECTION 

A principal advantage of treating the degenerate case, where the states involved in the 
matrix elements have the same energy, is the simplicity of the formalism. In the more general 
case, when the energies are different, the radiative correction to the weak interaction matrix 
element has to involve the laser field photon that drives the transition, otherwise one would 
not be dealing with a gauge-invariant amplitude. In the degenerate case we can restrict 
our attention to the gauge- invariant subset of diagrams shown in Fig. 1, which involve the 
vertex (Fig. lb) and wave function (Figs, la, lc) corrections. While the treatment of these 
diagrams is straightforward for scattering processes, more care is required when bound states 
are involved. As mentioned in the introduction, the similar problem of radiative correction 
to hyperfme splitting has already been treated in the literature [4, 5] , but in the present case 
the initial and final states are different, and the formalism requires some modifications. 

The bound state wave functions ip v and ip w are solutions of the Dirac equation in the field 
of a point nucleus. Therefore they can be interpreted as residues at poles of Dirac-Coulomb 
propagators as a function of energy E = p°, 

W^Wjz^V"^- (8) 

When radiative corrections are involved, the Dirac-Coulomb propagator is corrected by the 
electron self-interaction £ 

^-m-7°V-£(£)l r )" (9) 
The new position of the pole and corresponding residues are 

Ef = E^ + {4,\Z(E^) (10) 



~YdE 



+ ^(^)£(^M + 



E=E,i 



imm), (ii) 



where by S' F one denotes a reduced Coulomb-Dirac propagator, namely the propagator with 
the ^-state excluded. With the help of the above equations, we now present the one-loop 
radiative corrections to Q . They consist of the vertex correction Q v , the left and right 
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wave function corrections Qsl + Qsr, which include as well the derivative terms, associated 
with the last term in Eq. (11). In the Feynman gauge they are (e = E v — E w ) 

Q v = -Ama J fc2 i - (w\e ik -%S F (e - k )-f -f 5 p N S F (e - ko^e'^v), (12) 
f d n k 1 

Q SL = -Ama J (w\iol5pNS' F (e)e lk ^^S F (e - ko^e- 1 ^} 

+2««g |^^^( W |e l %^(6-A; )7o^(e-A; )7^- lfc% >) ) (13) 

r d n k 1 

Q SR = -Ama J (w\e lk ^^S F (e - k )Ye- ik ^S' F (e) lol5 p N \v) 

r d n k 1 

+2maQ j ^ ( v \ e ^S F (e - k ) lo S F (e - k )Ye- lk ^\v) . (14) 

There is still an ambiguity in the above formulas, related to the fact that at least one of the 
states is unstable with respect to radiative decay. This means that, for example, derivative 
terms, which have the interpretation of bound state wave function renormalization, acquire 
a small imaginary part. We think that this imaginary term may have a small effect on the 
weak matrix element. Nevertheless, in our treatment we completely ignore this imaginary 
part for simplicity. To include it properly would require a more detailed treatment of the 
excitation and decay process. Before the numerical integration, we present in the next 
section the analytic calculation of the first two terms in the Za expansion. 

IV. Za EXPANSION 

In the Za expansion one performs a simplification, similar to that used for the Lamb 
shift, which leads to an exact expression for the expansion terms. Specifically, the first two 
terms are given by the on-mass-shell scattering amplitude, which because it involves the 
weak charge of the nucleus, is dominated by the large momentum region, with characteristic 
momenta of the order of the electron mass, and to smaller extent of the order of the inverse of 
nuclear size. The small momentum region contributes at order O(Za) 2 and will be included 
in the numerical treatment. We aim here to confirm the previously obtained result [12] 
shown in Eq. 3, which will be used later to test the numerical accuracy of the nonperturbative 
treatment. In this section we do not pull out a factor a/ir from R(Za). 
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The relative correction to order a is determined by considering the radiative correction 
to the 7 M 7 5 vertex, 

tw \ = a f rf4 g W(P2,Pi) , . 

\P2,Pi) A7r J i7r 2^ q _ k y_ m 2 + ie ^ q 2_ m 2 + ie ^ q _ p2 y_ X 2 + ie y ^> 

where 

N»(p2,Pi) =l a (4- # + m) 7 VU+™ha, (16) 
and k = p 2 — pi- The most general form of T M in momentum space is 

r^(p 2 , Pl ) = F 1 (k 2 )-f'f + ^ 2 (A; 2 )-7 5 . (17) 

The form factors Fi(k 2 ) and F 2 (k 2 ) are calculated following the same steps as in the case of 
the electromagnetic vertex. Introducing Feynman parameters and taking into account the 
mass-shell condition, one obtains in the limit of zero momentum transfer 

fi(0) = ~ (18) 

For a static nucleus, only -F\(0) contributes to the relative correction to first order in a. The 
relative correction to the PNC amplitude is 

= uMr%M 

u{p, o r j7 u 7°'u(p, a - J 



which can be transformed into 



Tr 
Tr 



= lro^ +m) (i + ^ 

11 7°7 5 ik(^ + m )( 1 + ^ 5 ) 



where a M = (ao,a) with a • p = is the polarization four vector of the electron. We then 
recover the well-known [11] lowest-order correction 

R(Za) = (22) 

The leading binding correction can be derived from the forward scattering amplitude, 

which involves an additional Coulomb exchange. It consists of the 4 diagrams presented in 

Figs. 2a - 2d, which we evaluate using Yennie gauge. This gauge has the useful property 

that each diagram is infrared finite as the photon mass A is taken to 0. The contribution 

from Fig. 2a to the ratio R(Za) can be written as 

_Za^ r _d^l_ r d*q N^k) 

1_ a J {2nfk 2 J n 2 i [q 2 + ie] 2 [(p + q) 2 - m 2 + ie} 2 [(p + k + q) 2 - m 2 + ie}' 1 ' 
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where 



Ni(q,k) = (g^q 2 + 2g^g,)Tr[ 7 ^(^+^+m)7 (^+^+^ + m) 7 7 5 
x {$ + d + m ) 7 "_L(^ + m)(l + ^ 7 5 ) 

For Figs. 2b and 2c, we have 

Za 2 f d 3 k 1 fd 4 q 1 



where 



Za" r d°k 1 r d 



x 



2 i [g 2 + ie] 2 

Ni(q,h) 



[(p + q) 2 — m 2 + ie] [{p + k + q) 2 — m 2 + ie] [(p + k) 2 — m 2 + ie] ' 



(24) 



(25) 



N 2 (q,k) = (^g 2 + 2g^)Tr[ 7 u (^+^ + m)7^(^+^ + ^+m)7 7 5 
x + ^ + m)7^-^(^ + m) (1 + ^7 5 ) 



(26) 



W 3 (?,fc) = (^g 2 + 2g M g,)Tr[7^(^+^+m)7 (^+^ + ^+m) 7 " 



x # + m)7°7 5 — (jf + m)(l + ^7 5 ) 

^ 1 1 L 



Finally, for Fig. 2d, one has 

Za 2 /■ rf 3 A; 1 fd A q 



Ra — 



N 4 (q,k) 



(27) 



a J (2tt) 3 k 2 J n 2 i [q 2 + ie] 2 [(p + k) 2 — m 2 + ie] 2 [(p + k + q) 2 — m 2 + ie] 



(28) 



Ni{q,k) = (g fluq 2 + 2q fl q u )TT[ 1 {p'+^ + m)Y{^+4+m)Y{^+^ + m) 1 
x 7 5 ^(^ + m)(l + ^7% 

where k = (0, fc). Each contribution from Figs. 2a, 2b, 2c and 2d is written as 

Za 2 r d 3 k F t (k 2 ) 



(29) 



a J (2tt) 3 k 2 ■ (30) 

The calculations are considerably simplified if one determines only the imaginary part of the 
functions Fi(k 2 ). These are analytic functions with a branch cut for k 2 > 0. The real part 
of Fi(k 2 ) is then obtained by means of Cauchy's theorem, 



yj 2mJ M 2 -k 2 it J M 2 -k 2 



(31) 
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where k 2 < 0. Substituting this expression into Eq. (30) and integrating over k yields 

R * = i^-J { Q dM%F t (M 2 )}. (32) 

In order to calculate S[-Fi], a procedure in Mathematica is written which facilitates the 
evaluation of the trace in Eq. (24) and the the integrals in Eq. (32). Each contribution is 
doubled due to the permutation of photon and boson lines. Setting m—1 and picking the 
terms linear in p, we obtain 

A^rw^M f7 32 2 1 / 16 1 M\ 

g^(M )] = 2a o7 r(--^ + 3(i + M2)2 - TT ^ + 2^--- y j 



x 



2 



arctan(M) - arccos(^(M - 2)j - 2 [l - ^) Jl - A 6{M - 2) J, (33) 



where is the step function with 9(x) = for rr < and = 1 for x > 0. The above 
expression can be analytically integrated. Hence we have 

R(Za) = _« + ^ = _- + g. J ^ (Uf ^ a[F(] = _^_(^ +21n2 ) Za2 , ,34) 
in agreement with Ref. [12]. We now turn to the numerical calculation. 



V. NUMERICAL APPROACH 

In order to make contact with the notation used in Ref. [5], we note that the two terms 
Qsl and Qsr in Eqs. (13) and (14) are associated with what are called "side-left" (SL) 
and "side-right" (SR) diagrams in that work, which notation we will follow in this section. 
In addition, the SL and SR diagrams have contributions called "derivative terms" . We will 
refer in this section to the Gell-Mann Low formalism used in Ref. [5] in a rederivation of Eqs. 
(13) and (14): the adiabatic damping factor e used in that formalism can be distinguished 
from the factor used in dimensional regularization, n = 4 — e, by context. In the numerical 
evaluation, each diagram breaks into several pieces, which we define as 

Q = E Qv* + E Qsu + E Qsm- (35) 

i=l i=l i=l 

We now treat the vertex and side diagrams in turn. 
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A. Vertex diagram 



The vertex diagram Q v , shown in Fig. lb, was given in Eq. (12). The ultraviolet divergent 
part of the diagram can be isolated by replacing Sp with So, where So is a free propagator. If 
this replacement is made, we get the contribution Q V1 which is most conveniently evaluated 
in momentum space 

d 3 p 2 d 3 p Y r d n k 1 -r 1 „, , 1 



Qvi = —4:iria 



(2tt) 3 (2tt) 3 J (2vr)« k 2 + iS 



— ft — m 



V(q)- 



ii — ft — m 



(36) 



After Feynman parameterization the d n k integration can be carried out with the result 

a / C 



Qvi — 



2vr V e 



i )Qo 



/ pdp dx 
Jo Jo J 

d 3 p 2 d 3 p! 



11 j f d 3 p 2 d 3 p Y T Ay 



(2tt) 3 (2tt) j 



a f 1 f 1 , f d 3 p 2 d 3 p 
-A-Jo pdp L dX lj2nfJ^f 



x 



(37) 



Here 



C = (47r) e / 2 r(l + e/2), 
<5 M = pxp lfl + p(l - x)p 2fl , 
Ay = px(m 2 - p\) + p(l - x)(m 2 - p 2 ) + Q 2 , 
Q = \P2-Pi\, 

and the Fourier transform of the weak Hamiltonian in the case of a uniform charge distri- 
bution is 



V(q) 



sin(gi? ) — qRoCOs(qRo) 



7o75- 



(38) 



8n 3 (qR ) 3 

The first two terms in the right-hand-side of Eq. (37) are divergent and will be held for later 
cancellation with the "derivative terms" from the SL and SR calculation. The remaining 
finite parts of Q V1 are tabulated in the second column of Table II. 

The difference of Qy and Q V1 is ultraviolet finite, and is evaluated in coordinate space. 
The k integral is treated by carrying out a Wick rotation, k — > iuj, which leads to 

r°° duo f d 3 k e^-t) 



Q V2 = Sir a $tjd 3 xj d 3 y J d 3 z j™ J- 



x 



(2vr) 3 ^ + k 2 

i) w {x)-i^S F (x, y, e w - i^hol5pN(y)S F (y, z; e v - iu^i/j^z) 



- il> w (x)'y li So(x, y, e w - iuj)iol<opN{y)S Q (y, z; e v - iuj)^ip v (z) 



(39) 
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A singularity associated with the parts of the bound propagators which include w or v is 
regularized by evaluating the expression with e w<v — > e W}V (l — A) in the electron propagators. 
The result behaves as In A. While it is possible to explicitly cancel this dependence with 
similar terms from the side diagrams, we choose here to simply work with a specific, small 
value of A = 1CT 5 . We note that different choices of A will lead to slightly different results 
for Qv2, but when combined with the side diagrams discussed below, the sums are essentially 
the same as long as the values of A are reasonably small. Results for Q V 2 with A = 1CT 5 
are given in the third column of Table 11. 

The Wick rotation mentioned above passes bound state poles which must be accounted 
for. They are treated by rewriting Qy by treating the propagators as a spectral representa- 
tion, carrying out the d 3 k integration analytically, and defining 

/ e iVE 2 +i8\x-y\ 
d 3 xd 3 y — — — — — ^(x^^Mx) ^(yh^Mv), (40) 
\x y 



which allows us to write 



. r dko t—^ 9wnmv{k<i)Qmn ( AA\ 

i ~2tt ^ [e w (l — A) — k — e m (l - i5)} [e„(l - A) - k - e n (l - iS)] ' 1 ' 
The choice we have made in regularizing leads to only the ground state lsi/2, denoted as a, 
being encircled when k — > icu, so 

-e a )Q ma i . n \ 

W = L + • l 42 J 

an t a fc n ma t a t m 

The sum over a ranges only over the two magnetic quantum numbers of the state. This 
contribution is tabulated in the fourth column of Table II. The part of the summation in 
which the denominator would vanish corresponds to a double pole, but does not contribute 
because Q aa vanishes. However, it should be noted that double poles will in general con- 
tribute, and in fact would be present in the present calculation were we to use a negative 
value of A which would introduce additional pole terms from the 2si/ 2 and 2p 1 / 2 states. 

B. Side diagrams 

It is convenient for the discussion of the side diagrams to introduce the matrix element 
of the self-energy operator between two arbitrary states m and n, 

/f d n k ^-(x-y) 
d 3 xd 3 y J j—^ -^-^ $ m {x)^S F (x, y-E- k^Mv)- (43) 
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A self-mass counterterm is understood to be included in the above. The self-energy of a 
valence state is then £ TO (e„), and can be evaluated as described in Ref. [16]. 

Using a spectral decomposition of the intermediate propagator, the S'-matrix for SL is 

Qsl = -*A 3 £ / / d ^ Q^UE 2 ) _ _ ^ _ ko)A ^ + ^ _ 

m J Ztc J 2tc E\ — e m (l — 10) 

(44) 



and for SR 



Qsr = -*A 3 E / ^ / ^ ^ E ^ Qm ] ' A(E 2 - e v )A(E 1 - E 2 + k Q )A(E 1 + k - e w ), 
m J 2ir J 2-k E 2 — e m [l — id) 

(45) 

with 

A{E) = wh- (46) 

Here A is a factor associated with the Gell-Mann-Low formalism [17] that is to be differen- 
tiated and set to unity: in addition, a factor ie/2 must be multiplied into the S'-matrix to 
obtain the off-diagonal energy. If the restriction is made that m ^ w,v, it is straightforward 
to show that two "perturbed orbital" (PO) contributions to the matrix element result which 
are given by 

Qsli = E QwmJ] ™ {ev) ee £,„(e„) (47) 

and 

vsiji — 2^ ~ = ^wwK^w)- (4° J 

This is equivalent to the forms given for Qsl and Qsr in Eqs. (13) and (14). We note that 
it is not necessary to explicitly make the restrictions m ^ w in Qsli and m ^ v in Qsri 
because Q ww = Q vv = 0. The PO terminology arises from the fact that the m summation 
can be carried out before evaluating the self-energy, and one then needs only to do a self- 
energy calculation with one of the external wavefunctions replaced with a perturbed orbital. 
The PO terms are tabulated in the fifth and ninth columns of Table II. 

The cases m = v and m = w are more subtle, as they contribute terms of order 1/e to the 
off-diagonal energy. This divergence cancels, but a finite contribution coming from Taylor 
expanding T,(E) remains, and contributes |Q ^w(et,) + |QoS' U)W (e w ), or more explicitly 

d n k e^-^ 



Qf[ = 2maQ J d 3 yd 3 zd 3 w J 



(2tt)™ k 2 + iS 

x4>v(yhnS F (y, w; e v - k )>y S F (w, z; e v - k )>y^v(z) 
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2 ^ J 2-k (e v - k - E m y 



and 



Q% = 2iiraQ J d 3 yd 3 zd 3 w J 



d n k e ik - ( y- 2) 



-2 Q »\?]^(e„,-k n -e r (50) 



(2tt)™ k 2 + iS 
x^w1^S F (y, w; e w - /c )7o5'f(w, z\ e w - k )^ w (z) 

i ^ \ ^ f dk() Qwmmwiko) 

These correspond to the derivative term mentioned in the preceding section. The analysis 
of this term parallels closely the treatment of the vertex: first the bound propagators are 
replaced with free propagators, which gives 

/dp f d n k — 1 1 

(2^7 J^MP^JI^^^j^^rMp) (51) 

and 

/d^p f d n k — 1 1 

J^fl J2nf ^^J^J^^J^J^ 1 "^ (52) 
Feynman parameterizing and carrying out the d n k integration gives 



2 



Qsl2 = -^(7 - i)Qo + £q f Q pd P I -^Mp) 1o Mp)^ As 

+ ^ Qo Jo pdp JW 



a Q fjdp J J^Pfaffi^Mi - p ) + m] 7 o[^(l - P) + m] 7 ^(p)}^- (53) 



and 



Q SR 2 = -^Qo + ^Qoj\dpf^jMloMp)^ 

/ (2^3 {^^t^ 1 ~ p)+ m ^oW ~P)+ m] 7 ^(p)}^-(54) 
where A 5 = p 2 p 2 + p(m 2 — p 2 ). The first two terms in the right-hand-side of the above 
equations for Qsl2 and QsR2 are divergent but cancel with the corresponding terms of Qvi 
in Eq. (37). The remaining, finite terms are presented in the sixth and tenth columns of 
Table II. 

The difference of the side diagrams evaluated with bound propagators and free propaga- 
tors is again ultraviolet finite, and after Wick rotation one has 

'duj f d 3 k e lk i s ~z) 



(27r) 3 u 2 + k 2 

$v{x)1vlS f {x, y, e v - iu)~f G S F (y, z\ e v - iu)^i) v {z) 
- ^ v (x)^S (x, y, e v - iu)^ S (y, z\ e v - iu)^ip v (z) 



(55) 
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and 

x ^ w (x)lf^S F (x, y, e w - iu)^ Q S F {y, z\ e w - iu)Y^w{z 
- rl> w (x)'y l J5o(x , y, e w - iuj)^S {y, z; e w - iu>)^ip w (z) 



(56) 



The same regularization of the valence energy, — > e„(l — A) used in the vertex is required, 
and again we simply use A = 1CT 5 and present the results in the seventh and eleventh 
columns of Table II. 

Finally, the Wick rotation passes a double pole when m = a, with a being the lsi/2 
ground state, leading to the derivative terms 

Qsla = \Qo ]T g' vaav (e v - e ) (57) 

and 

Qsm = lQoY,9' waaw ( e w ~ e a ) (58) 

which are tabulated in the eighth and twelfth column of Table II. This completes the 
calculation and the sums of vertex and side diagram contributions give the values of the 
exact evaluation of the function R(Za) which are tabulated in the last column of Table II. 



VI. DISCUSSION 

A number of numerical issues arise in the calculation that we note here. In some parts 
the use of a uniform distribution, with its step function behavior, caused loss of accuracy. In 
those cases a fermi distribution was used: while this leads to small changes in Q , the effect 
on R(Za) is negligible. More serious is the difficulty of controlling numerical instabilities 
at low Z, which led to our choosing the lowest Z to be 10. A graph of the numerical value 
of R(Za), along with results from the two leading terms given in Eq. 3, is shown in Fig. 3. 
The accuracy of the calculation at low Z is sufficient to allow a fit that determines 

R(Za) Bt = ~-0M5(2)Z. (59) 
which is in agreement with the numerical value of Eq. (3) 

R(Za) = ~\- 0.045154Z. (60) 
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This agreement provides a check on the rather complex numerical calculation. The ad- 
vantage of the numerical approach is of course the fact that it does not assume Za to be 
a small parameter, and thus can be used for high Z. In the particularly interesting case of 
Z = 55, we see that the perturbative formula happens to be -2.983, as compared with the 
exact result -4.007. 

However, at this point we make no claims about the applicability of the present calcu- 
lation to the case of PNC transitions in neutral cesium. The actual process studied in the 
experiment that measures PNC [6] involves a double perturbation, where not only the weak 
Hamiltonian but also an external laser photon field act to either first transform the 6si/2 
electron into a state with the opposite parity, followed by an allowed dipole transition to a 
7si/2 electron, or vice- versa. To extend our calculation to the actual experiment requires 
the following steps. 

The first step is replacing the Coulomb wave functions used here with realistic wave 
functions for neutral cesium. The technology to carry out radiative corrections in neutral 
atoms has only recently been put into place. It is now possible, using a local potential that 
incorporates screening, to carry out accurate self-energy [13, 14] and radiative correction to 
hfs [15] calculations. The second step is to incorporate the laser photon. This is a more 
complicated task, since the set of diagrams shown in Fig. 4 must be evaluated. We note that 
while Fig. 4a corresponds to the vertex correction considered here, Fig. 4c corresponds to a 
radiative correction to the electromagnetic vertex, and Fig. 4e to a new radiative correction 
specific to the experiment. We expect that the new contributions will affect R(Za) in order 
(Za) 2 , but until they are explicitly evaluated, their importance for cesium PNC is unknown. 

The principal results of this paper are then as follows. Firstly, an independent scattering 
calculation of the leading binding correction in the function R(Za) has been presented, 
which confirms the calculation of Ref. [12]. Secondly, it has been shown that numerical 
methods that allow the evaluation of R(Za) to all orders in Za for the case of gauge- 
invariant 2si/ 2 — 2p!/ 2 transitions in hydrogenlike ions can be applied. The behavior of the 
function shows that large binding corrections are present. If these corrections behave the 
same in the realistic cesium case, and the extra corrections of order (Za) 2 mentioned above 
are small, the apparent 2a discrepancy noted in Ref. [8] is reduced, but before a complete 
calculation is completed the theoretical status of PNC in cesium should be regarded as 
unresolved. 
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I I I 

XXX 
(a) (b) (c) 



FIG. 1: Feynman diagrams for the self-energy corrections to parity nonconservation. The dashed 
line terminated with a cross indicates an interaction with the nucleus through the exchange of a Z 
boson. 




(a) (b) (c) (d) 

FIG. 2: Feynman diagrams for the leading binding corrections to the forward scattering amplitude. 
The dashed and wavy lines terminated with a cross represent PNC and Coulomb interactions with 
the nucleus, respectively. 




Z 

FIG. 3: Comparison of exact calculation (solid line) with first two terms of perturbative expansion 
for R(Za) (dashed line). 
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FIG. 4: Feynman diagrams for the radiative correction to electron excitation by a laser photon, 
indicated by the wavy line terminated with a triangle, in the presence of interaction with the 
nucleus through exchange of a Z boson, indicated by the dashed line terminated with a cross. 
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TABLE I: Nuclear parameters c and i?o and lowest-order PNC matrix element Qq: units of fermis 
for c and Rq and 1/ajj for Qq. Square brackets indicate power of 10. 



7 
/ . 


c 


ito 


dr. 
WO 


1 n 


z.yooy 


o.ooy 


l.olo[UJ 


10 


o.Z ( 0Z 


/l 1 97 
4.1Z 1 


f .Uoo[U| 


9D 
ZU 


o. i loo 


A A%7 
4.4o 1 


9 ^8R[1 1 
Z.OOO[lJ 


9^ 
ZO 


4.U f UO 


A 7SQ 
4. f OO 


D.ODD[lJ 


ou 


A AA <^A 


0. 1UD 


1.400 |ZJ 


/in 

4U 


a oi 1 c; 
4.y i io 


0.01D 


0.yoo|Zj 


ou 


o.4o»o 


r ni n 

D.U1U 


Z.UlU|oJ 


55 


5.6748 


6.206 


3.539[3l 

L J 


60 


5.8270 


6.345 


6.136[3] 


70 


6.2771 


6.761 


1.786 [4] 


80 


6.6069 


7.068 


5.184[4] 


90 


6.9264 


7.368 


1.542 [5] 


100 


7.1717 


7.599 


4.886[5] 



19 



TABLE II: Breakdown of Contributions to R{Za). 



z 


Qvi 


QV2 


QV3 


Qsli 


QSL2 


QSL3 


QsLi 


Qsm 


QSR2 


Qsm 


Qsm 


R(Za) 


10 


-2.500 


4.260 


-11.768 


-0.272 


2.369 


2.293 


0.002 


-0.045 


2.269 


2.388 


0.006 


-0.998 


15 


-2.009 


-0.253 


-7.853 


-0.415 


1.972 


2.689 


0.003 


-0.068 


1.873 


2.782 


0.009 


-1.270 


20 


-1.729 


-2.600 


-5.898 


-0.555 


1.694 


2.967 


0.004 


-0.097 


1.595 


3.058 


0.012 


-1.549 


25 


-1.559 


-4.064 


-4.727 


-0.696 


1.482 


3.180 


0.005 


-0.128 


1.384 


3.269 


0.014 


-1.840 


30 


-1.458 


-5.067 


-3.948 


-0.839 


1.311 


3.350 


0.006 


-0.164 


1.214 


3.439 


0.017 


-2.139 


40 


-1.377 


-6.398 


-2.980 


-1.139 


1.048 


3.620 


0.008 


-0.279 


0.952 


3.703 


0.022 


-2.820 


50 


-1.386 


-7.267 


-2.406 


-1.461 


0.851 


3.824 


0.011 


-0.434 


0.755 


3.904 


0.026 


-3.583 


55 


-1.411 


-7.609 


-2.201 


-1.637 


0.768 


3.911 


0.012 


-0.530 


0.673 


3.989 


0.028 


-4.007 


60 


-1.446 


-7.912 


-2.033 


-1.821 


0.694 


3.990 


0.013 


-0.644 


0.600 


4.067 


0.030 


-4.462 


70 


-1.525 


-8.466 


-1.778 


-2.229 


0.566 


4.132 


0.016 


-0.922 


0.472 


4.207 


0.033 


-5.494 


80 


-1.617 


-8.922 


-1.602 


-2.707 


0.458 


4.253 


0.019 


-1.300 


0.365 


4.331 


0.035 


-6.687 


90 


-1.708 


-9.376 


-1.489 


-3.288 


0.364 


4.366 


0.023 


-1.822 


0.273 


4.448 


0.035 


-8.174 


100 


-1.798 


-9.831 


-1.433 


-4.020 


0.282 


4.457 


0.028 


-2.570 


0.192 


4.545 


0.033 


-10.115 



20 



